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Abstract 

The present experimental status in the search for neutrinoless double 
beta decay is reviewed, with emphasis on the first indication for neutri- 
noless double beta decay found in the HEIDELBERG-MOSCOW exper- 
iment, giving first evidence for lepton number violation and a Majorana 
nature of the neutrinos. Future perspectives of the field are briefly out- 
lined. 



1 Introduction 

The neutrino oscillation interpretation of the atmospheric and solar neutrino 
data, deliver a strong indication for a non-vanishing neutrino mass. While such 
kind of experiments yields information on the difference of squared neutrino 
mass eigenvalues and on mixing angles (for the present status see, e.g. [54, 55, 
56, 57]), the absolute scale of the neutrino mass is still unknown. Information 
from double beta decay experiments is indispensable to solve these questions 
[5, 7, 8]. Another even more important problem is that of the fundamental 
character of the neutrino, whether it is a Dirac or a Majorana particle [9, 10]. 
Neutrinoless double beta decay could answer also this question. Perhaps the 
main question, which can be investigated by double beta decay with highest 
sensitivity, is that of lepton number conservation or non-conservation. 

Double beta decay, the rarest known nuclear decay process, can occur in 
different modes: 

2v/3/3 - decay : A(Z, N) — > A(Z+2, N—2) + 2er + 2v e (1) 
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Oz//3/3 - decay : A(Z, N) A(Z+2, N-2) + 2e~ (2) 

0u(2) X P(3 - decay : A(Z, TV) — > A(Z+2, N-2) + 2eT + (2)x (3) 



While the two-neutrino mode (1) is allowed by the Standard Model of particle 
physics, the neutrinoless mode (Qvj3(3 ) (2) requires violation of lepton number 
(AL=2). This mode is possible only, if the neutrino is a Majorana particle, 
i.e. the neutrino is its own antiparticle (E. Majorana [9], G. Racah [10], for 
subsequent works we refer to [12, 13, 14], for some reviews see [26, 27, 28, 29, 8]). 
First calculations of Qv{3j3 decay based on the Majorana theory have been done 
by W.H. Furry [11]. The most general Lorentz-invariant parametrization of 
neutrinoless double beta decay is shown in Fig. 1. 



Figure 1: Feynman graphs of the general double beta decay with long range 
(a-c) and short range (d) parts (see [16]). 

The usually used assumption is that the first term (i.e. the Majorana mass 
mechanism) dominates the decay process. However, as can be seen from Fig. 
1, and as discussed elsewhere (see, e.g. [8, 32, 33, 47]) neutrinoless double 
beta decay can not only probe a Majorana neutrino mass, but various new 
physics scenarios beyond the Standard Model, such as R-parity violating super- 
symmetric models, R-parity conserving SUSY models, leptoquarks, violation of 
Lorentz-invariance, and compositeness (for a review see [8, 32, 33]). Any the- 
ory containing lepton number violating interactions can in principle lead to this 
process allowing to obtain information on the specific underlying theory. It has 
been pointed out already in 1982, however, that independently of the mecha- 
nism of neutrinoless double decay, the occurence of the latter implies a non-zero 
neutrino mass and vice versa [15]. This theorem has been extended to super- 
symmetry. It has been shown [17] that if the neutrino has a finite Majorana 
mass, then the sneutrino necessarily has a (B-L) violating 'Majorana' mass, too, 
independent of the mechanism of mass generation. 




u 
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The experimental signature of the neutrinoless mode is a peak at the Q-value 
of the decay. 

Restricting to the Majorana mass mechanism, a measured half-life allows to 
deduce information on the effective Majorana neutrino mass (m) , which is a 
superposition of neutrino mass eigenstates [26, 27]: 

PT/ 2 (o+ - o+r 1 = c m J-^- + c nv { v ) 2 + c xx (\f + c mv { v ) { -^ 

+C mX (X)^- + C nX ( V )(X), (1) 

7Yl e 

( m ) = | m a)| + e ^| m (2)| + e ^3| m (3)| ; (2) 

where mi 1 } = \rnie | exp (i4>i) (i = 1,2,3) are the contributions to (m) from in- 
dividual mass eigenstates, with <fri denoting relative Majorana phases connected 
with CP violation, and C mm ,C w ,... denote nuclear matrix elements squared, 
which can be calculated, (see, e.g. [39], for a review and some recent discussions 
see e.g. [40, 8, 27, 28, 43, 45, 46, 47]). Ignoring contributions from right-handed 
weak currents on the right-hand side of eq.(l), only the first term remains. 
The effective mass is closely related to the parameters of neutrino oscillation 



experiments, as can be seen from the following expressions 

|m«| = |C/ e i| 2 mi, (3) 

\m%\ = \U e2 \*yfAml 1+ ml (4) 

\m%\ = |[/ e3 |yAm| 2 + Am2 1+ m?, (5) 

Here, U e i are entries of the neutrino mixing matrix, and Amfj = \mf — m||, 



with nii denoting neutrino mass eigenstates. U e i and Am 2 can be determined 
from neutrino oscillation experiments. 

The importance of (m) for solving the problem of the structure of the neu- 
trino mixing matrix and in particular to fix the absolute scale of the neutrino 
mass spectrum which cannot be fixed by v - oscillation experiments alone, has 
been discussed in detail in e.g. [5, 7]. 

Double beta experiments to date gave only upper limits for the effective 
mass. The most sensitive limits [19, 20, 21] were already of striking importance 
for neutrino physics, excluding for example, in hot dark matter models, the 
small mixing angle (SMA) MSW solution of the solar neutrino problem [48, 49, 
50, 51, 52, 53, 7, 8] in degenerate neutrino mass scenarios. 

The HEIDELBERG-MOSCOW double beta decay experiment in the Gran 
Sasso Underground Laboratory [4, 18, 23, 25, 32, 24, 8, 21] searches for double 
beta decay of 76 Ge — > 76 Se + 2 e~ + {2v) since 1990. It is the most sensitive 
double beta experiment since nine years now. The experiment operates five 
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enriched (to 86%) high-purity 76 Ge detectors, with a total mass of 11.5 kg, the 
active mass of 10.96 kg being equivalent to a source strength of 125.5 mol 76 Ge 
nuclei. 

We have performed recently [1, 2, 3] a refined analysis of the data obtained 
during the period August 1990 - May 2000 [21]. This analysis concentrates 
on the neutrinoless decay mode which is the one relevant for particle physics 
(see, e.g. [8]), and reports first evidence for the neutrinoless decay mode (see 
also [110, 111, 112, 113, 114]). First reactions have been published [5, 58, 59, 
60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80] 
etc.). For the results concerning 2v/3(3 decay and Majoron-accompanied decay 
we refer to [21, 116]. We shall describe here the results leading to the evidence 
for 0v(3(3 decay. The results will then be put into the context of other present 
f}(3 activities. 

2 Experimental Set-up and Results 

A detailed description of the HEIDELBERG-MOSCOW experiment has been 
given recently in [18], and in [21]. Therefore only some important features will 
be given here. We start with some general notes. 
1. Since the sensitivity for the Qv(3(3 half-life is 



(and ~ (m„)), with a denoting the degree of enrichment, e the efficiency 

of the detector for detection of a double beta event, M the detector (source) 
mass, AE the energy resolution, B the background and t the measuring time, 
the sensitivity of our 11kg of enriched 76 Ge experiment corresponds to that 
of an at last 1.2 ton natural Ge experiment. After enrichment, energy resolu- 
tion, background and source strength are the most important features of a f}(3 
experiment. 

2. The high energy resolution of the Ge detectors of 0.2% or better, as- 
sures that there is no background for a Qv(5(3 line from the two-neutrino double 
beta decay in this experiment, in contrast to most other present experimental 
approaches, where limited energy resolution is a severe drawback. 

3. The efficiency of Ge detectors for detection of Qv(3j3 decay events is close 
to 100% (95%, see [81]). 

4. The source strength in this experiment of 11kg is the largest source 
strength ever operated in a double beta decay experiment. 

5. The background reached in this experiment is, with 0. 17 events/kg ykeV 
in the 0v(3[3 decay region (around 2000-2080 keV), the lowest limit ever obtained 
in such type of experiment. 

6. The statistics collected in this experiment during 10 years of stable run- 
ning is the largest ever collected in a double beta decay experiment. 




(6) 
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Figure 2: The /3/3-laboratory of the HEIDELBERG-MOSCOW experiment in 
the Gran Sasso and its location between halls A and B, and four of the enriched 
detectors during installation (from up to bottom). 
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Total 


Active 


Enrichment 


FWHM 1996 


FWHM 2000 


Detector 


Mass 


Mass 


in Yb Ge 


at 1332 keV 


at 1332 keV 


Number 


[kg] 


[kg] 


[%] 


[keV] 


[keV] 


No. 1 


0.980 


0.920 


85.9 ± 1.3 


2.22 ± 0.02 


2.42 ± 0.22 


No. 2 


2.906 


2.758 


86.6 ± 2.5 


2.43 ± 0.03 


3.10 ± 0.14 


No. 3 


2.446 


2.324 


88.3 ± 2.6 


2.71 ± 0.03 


2.51 ± 0.16 


No. 4 


2.400 


2.295 


86.3 ± 1.3 


2.14 ± 0.04 


3.49 ± 0.24 


No. 5 


2.781 


2.666 


85.6 ± 1.3 


2.55 ± 0.05 


2.45 ± 0.11 



Table 1: Technical parameters of the five enriched detectors. (FWHM - full 
width at half maximum.) 

7. The Q value for ncutrinolcss double beta decay has been determined 
recently with very high precision to be Q ( 3 ( g=2039. 006(50) keV [82, 83]. 

We give now some experimental details. All detectors (whose technical pa- 
rameters arc given in Table 1 (see [18])), except detector No. 4, are operated 
in a common Pb shielding of 30 cm, which consists of an inner shielding of 10 
cm radiopure LC2-grade Pb followed by 20 cm of Bolidcn Pb. The whole setup 
is placed in an air-tight steel box and flushed with radiopure nitrogen in order 
to suppress the 222 Rn contamination of the air. The steel box is centered in- 
side a 10-cm boron-loaded polyethylene shielding to decrease the neutron flux 
from outside. An active anticoincidence shielding is placed on the top of the 
setup to reduce the effect of muons. Detector No. 4 is installed in a separate 
setup, which has an inner shielding of 27.5 cm electrolytical Cu, 20 cm lead, and 
boron-loaded polyethylene shielding below the steel box, but no muon shielding. 
Fig. 2 gives a view of the experimental setup. 

To check the stability of the experiment, a calibration with a 228 Th and a 
152 Eu+ 228 Th, and a 60 Co source is done weekly. High voltage of the detectors, 
temperature in the detector cave and the computer room, the nitrogen flow 
in the detector boxes, the muon anticoincidence signal, leakage current of the 
detectors, overall and individual trigger rates are monitored daily. The energy 
spectrum is taken in 8192 channels in the range from threshold up to about 3 
MeV, and in a parallel spectrum up to about 8 MeV. 

Because of the big peak-to-Compton ratio of the large detectors, external 7 
activities are relatively easily identified, since their Compton continuum is to a 
large extent shifted into the peaks. The background identified by the measured 
7 lines in the background spectrum consists of: 

1. primordial activities of the natural decay chains from 238 [7, 232 Th, and 

2. anthropogenic radio nuclides, like 137 Cs, 134 Cs, 125 Sb, 207 Bi 
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3. cosmogenic isotopes, produced by activation due to cosmic rays. 

The activity of these sources in the setup is measured directly and can 
be located due to the measured and simulated relative peak intensities of 
these nuclei. 

Hidden in the continuous background are the contributions of 

4. the bremsstrahlungs spectrum of 210 Bi (daughter of 210 Pb), 

5. elastic and inelastic neutron scattering, and 

6. direct muon- induced events. 

External a and (3 activities are shielded by the 0.7-mm inactive zone of the 
p-type Ge detectors on the outer layer of the crystal. The enormous radiopurity 
of HP-germanium is proven by the fact that the detectors No. 1, No. 2 and No. 
3 show no indication of any a peaks in the measured data. Therefore no contri- 
bution of the natural decay chains can be located inside the crystals. Detectors 
No. 4 and No. 5 seem to be slightly contaminated with 210 Pb on the level 
of few /xBq/kg, most likely surface contaminations at the inner contact. This 
contamination was identified by a measured a peak in the background spectrum 
at 5.305 MeV of the daughter 210 Po and the constant time development of the 
peak counting rate. There is no contribution to the background in the interest- 
ing evaluation areas of the experiment due to this activity. For further details 
about the experiment and background we refer to [18, 81] (see also Table 2). 

In the vicinity of the Q- value of the double beta decay of = 2039.006(50) 
keV, very weak lines at 2034.744 and 2042 kcV from the cosmogenic nuclide 
56 Co, and from 21i Bi ( 23S [/-decay chain) at 2010.7, 2016.7, 2021.8 and 2052. 9keV, 
may be expected. 

On the other hand, there are no background 7-lines at the position of an 
expected Ov(3(3 line, according to our Monte Carlo analysis of radioactive im- 
purities in the experimental setup [81] and according to the compilations in 
[84]. 

In total 55 possible 7-lines from various isotopes in the region between 2037 
and 2041 keV are known [84]. Only 5 of the isotopes responsible for them ( 102 Rh, 
146 Eu, 124 I, 124 Sb and 170 Lu) have half-lifcs larger than 1 day. However, some 
of the isotopes yielding lines in this energy range can in principle be produced 
by inelastic hadron reactions (induced by muons or neutrons). 

Therefore each of these 55 isotopes was checked for the existence of a 7- 
line from the isotope which has a high emission probability Ih- A search was 
made for this 7-line in the measured spectrum to obtain its intensity (Sh) or an 
upper limit for it. Then the adopted intensity So for a 7-line from the same 
isotope in the area around ^2039 keV can be calculated by using the emission 
probability Iq for the line at ~2039kcV. Different absorption for gammas of 
different energies are taken into account in a schematic way. 
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Detec- 
tor 


Life 
Time 


Date 


Shielding 


Background* 
[counts/ 
keV y kg] 


1 PSA 












boron- 


2000.-2100. 




Number 


[ days ] 


Start End 


Cu 


Pb 


poly. 




keV 




No.l 


387.6 


8/90-8/91 


yes 


0.56 


no 






1/92-8/92 






no 


No.2 


225.4 


9/91 - 8/92 


yes 


0.29 


no 


Common shielding for three detectors 


No.l 


382.8 


9/92 - 1/94 


yes 


0.22 


no 


No.2 


383.8 


9/92 - 1/94 


yes 


0.22 


no 


No.3 


382.8 


9/92 - 1/94 


yes 


0.21 


no 


No.l 


263.0 


2/94 - 11/94 




yes 


yes 




0.20 


no 


No.2 


257.2 


2/94 - 11/94 




yes 


yes 




0.14 


no 


No.3 


263.0 


2/94 - 11/94 




yes 


yes 




0.18 


no 


Full Setup 
Four detectors in common shielding, 


one 


detector separate 


No.l 


203.6 


12/94 - 8/95 




yes 


yes 




0.14 


no 


No.2 


203.6 


12/94 - 8/95 




yes 


yes 




0.17 


no 


No.3 


188.9 


12/94 - 8/95 




yes 


yes 




0.20 


no 


No.5 


48.0 


12/94 - 8/95 




yes 


yes 




0.23 


since 

2/95 


No.4 


147.6 


1/95 - 8/95 


yes 


0.43 


no 


No.l 


203.6 


11/95 - 05/00 




yes 


yes 




0.170 


no 


No.2 


203.6 


11/95 - 05/00 




yes 


yes 




0.122 


yes 


No.3 


188.9 


11/95 - 05/00 


yes yes 


0.152 


yes 


No.5 


48.0 


11/95 - 05/00 


yes yes 


0.159 


yes 


No.4 


147.6 


11/95 - 05/00 


yes 


0.188 


yes 



Table 2: Development of the experimental set-up and of the background num- 
bers in the different data acquisition periods for the enriched detectors of the 
HEIDELBERG-MOSCOW experiment. *) Without PSA method. 
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If the calculated intensity Si < 1 for the 7-line in the interesting area, this 
isotope can be safely excluded to contribute a significant part to the background. 

For example, the isotope 139 Xe possesses a 7-line at 2039.1 keV with an 
emission probability of 0.078%. This isotope also possesses a 7-line at 225.4 keV 
with an emission probability of 3.2%. The intensity of the line at 225.4 keV in 
our spectrum was measured to be <6.5 counts. This means that the 7-line at 
2039.1 keV has < 2^6.5 = 0.32 counts, and therefore can be excluded. 

Only eight isotopes could contribute a few counts according to the calculated 
limits, in the interesting area for the 0vf3(3 -decay area: 52m Fe, 93m Ru, 120 In, 
131 Ce, 170 Lu, 170m Ho, 174 Ta and 198 T1. Most of them have a half-life of a few 
seconds, only 170 Lu has a half-life of 2.01 days, and no one of them has a longer 
living mother isotope. To contribute to the background they must be produced 
with a constant rate, e.g. by inelastic neutron and/or muon reactions. Only 5 
isotopes can be produced in a reasonable way, by the reactions listed in Tabl.3. 



isotope 


52m-p e 


93m Ru 


120 In 


170m Ho 


198 rpj 


production 
reaction 


50 Cr(a, 2n 7 ) 


92 Mo(a, 3n) 


120 Sn(n,p), 
123 Sb(n,a) 


170 Er(n,p) 


197 Au(a, 3n 7 ) 



Table 3: Possible reactions (second line) which could produce the isotopes listed 
in the first line (see [88]). 

Except 120 Sb each of the target nuclides is stable. All reactions induced 
with a-particles can be excluded due to the very short interaction length of 
a-particles. Two possibilities remaining to explain possible events in the Qv(3f3 
-decay area would be: 

• 120 Sn(n,;p) 120 In: 

The cross section for this reaction is 2.5±1 mb for E n = 14.5 MeV [89]. 
Assuming a neutron flux of (0.4±0.4)xl0 -9 — K— for neutrons with an 
energy between 10-15 MeV as measured in the Gran Sasso [90] the rate 
of 120 In atoms produced per year is about 2 x 10~ 5 when there are 50 g 
of 120 Sn in the detector setup. Even when the cross-section is larger for 
lower energies, this can not contribute a significant number of counts to 
the background. 

• 170 Er(n,p) 170m Ho: 

The cross section for this reaction is about 1.13±lmb for E n — 14.8MeV 
[91]. Assuming again a neutron flux of (0.4±0.4)xl0~ 9 em ^ - for neutrons 
with an energy between 10-15 MeV the rate of 170 Er atoms produced per 
year is even less when assuming 50 g of 120 Er in the detector-setup. 
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In both cases it would be not understandable, how such large amounts of 
120 Sn or 170 Er could have come into the experimental setup. Concluding we do 
not find indications for any nuclides, that might produce 7-lines with an energy 
around 2039 keV in the experimental setup. 




energy [keV] 

Figure 3: Sum spectrum of enriched detectors Nr. 1,2,3,4,5 over the period Au- 
gust 1990 - May 2000 (54.9813 kg y) measured in the HEIDELBERG-MOSCOW 
experiment (binning 0.36 keV). The sources of the main identified background 
lines are noted. 

Fig. 3 shows the combined spectrum of the five enriched detectors obtained 
over the period August 1990 - May 2000, with a statistical significance of 54.981 
kgy (723.44 molyears). (Note that in Fig. 1 of [21]) only the spectrum of 
the first detector is shown, but normalized to 47.4 kgy [22]). The identified 
background lines give an indication of the stability of the electronics over a 
decade of measurements. The average rate (sum of all detectors) observed over 
the measuring time, has proven to be constant within statistical variations. 
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Figure 4: Sum spectrum of the 76 Ge detectors Nr. 1,2,3,4,5 over the period 
August 1990 to May 2000, (54.9813 kg y) in the energy interval 2000 - 2080 keV, 
around the Qpp value of double beta decay (Q/3/3 = 2039.006(50) keV) summed 
to 1 keV bins. The curve results from Bayesian inference in the way explained 
in sec.3. It corresponds to a half-life T^ 2 =(0.80 - 35.07) x 10 25 y (95% c.l.) 

on 4.5 r 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 11 
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Figure 5: Sum spectrum of single site events, measured with the detectors Nr. 
2,3,5 operated with pulse shape analysis in the period November 1995 to May 
2000 (28.053 kgy), summed to 1 keV bins. Only events identified as single site 
events (SSE) by all three pulse shape analysis methods [85, 86, 87] have been 
accepted. The curve results from Bayesian inference in the way explained in 
sec.3. When corrected for the efficiency of SSE identification (see text), this 
leads to the following value for the half-life: T^ 2 =(0.88 - 22.38) x 10 25 y (90% 



c.l.). 
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Fig. 4 shows the part of the spectrum shown in Fig. 3, in more detail around 
the Q-value of double beta decay. Fig. 5 shows the spectrum of single site 
events (SSE) obtained for detectors 2,3,5 in the period November 1995 - May 
2000, under the restriction that the signal simultaneously fulfills the criteria of 
all three methods of pulse shape analysis we have recently developed [85, 87] and 
with which we operate all detectors except detector 1 (significance 28.053 kg y) 
since 1995. 

Double beta events are single site events confined to a few mm region in the 
detector corresponding to the track length of the emitted electrons. In all three 
methods, mentioned above, the output of the charge-sensitive preamplifiers was 
differentiated with 10-20 ns sampled with 250 MHz and analysed off-line. The 
methods differ in the analysis of the measured pulse shapes. The first one relies 
on the broadness of the charge pulse maximum, the second and third one are 
based on neural networks. All three methods are 'calibrated' with known double 
escape (mainly SSE) and total absorption (mainly MSE) 7-lines [85, 86, 87, 81]. 
They allow to achieve about 80% detection efficiency for both interaction types. 

The expectation for a 0vf3/3 signal would be a line of single site events on 
some background of multiple site events but also single site events, the latter 
coming to a large extent from the continuum of the 2614keV 7-line from 2m Tl 
(see, e.g., the simulation in [85]). From simulation we expect that about 5% 
of the double beta single site events should be seen as MSE. This is caused by 
bremsstrahlung of the emitted electrons [81]. 

Installation of PSA has been performed in 1995 for the four large detectors. 
Detector Nr. 5 runs since February 1995, detectors 2,3,4 since November 1995 
with PSA. The measuring time with PSA from November 1995 until May 2000 
is 36.532 kg years, for detectors 2,3,5 it is 28.053 kg y. 

Fig. 6 shows typical SSE and MSE events from our spectrum. 




time [|ls] time [fts] 



Figure 6: Left: Shape of one candidate for Ov/3/3 decay classified as SSE by all 
three methods of pulse shape discrimination. Right: Shape of one candidate 
classified as MSE by all three methods. 

All the spectra are obtained after rejecting coincidence events between dif- 
ferent Ge detectors and events coincident with activation of the muon shield. 
The spectra, which are taken in bins of 0.36 keV, are shown in Figs. 4,5, Fig. 2 of 
[1] in 1 keV bins, which explains the broken number in the ordinate. We do the 
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analysis of the measured spectra with (Fig. 4) and without the data of detector 
4 (see Fig. 2 in [1], 46.502 kg y) since the latter does not have a muon shield 
and has the weakest energy resolution. The 0.36 keV bin spectra are used in all 
analyses described in this work. We ignore for each detector the first 200 days of 
operation, corresponding to about three half-lives of 56 Co (T 1 / 2 = 77.27 days), 
to allow for some decay of short-lived radioactive impurities. 

« 1000 



600 
400 



1430 1440 1450 1 460 1 470 1480 1400 

Energy [keV] 

Figure 7: Fit of one the lines (here the 1460.81 keV line from 40 K) in the sum 
spectrum of all five detectors, used for the calibration and the determination of 
the energy resolution (see Table 4). 



energy [keV] 
fit 


energy [keV] 
from [84] 


width [keV] 
fit 


width [keV] 
from calc. 


1460.81 ± 0.02 
1764.56 ± 0.05 
2103.31 ± 0.45 
2204.12 ± 0.14 
2447.73 ± 0.26 
2614.48 ± 0.07 


1460.81 
1764.49 
2103.53 
2204.19 
2447.86 
2614.53 


1.49 ± 0.01 
1.70 ± 0.05 
1.86 ± 0.35 
1.89 ± 0.13 
1.82 ± 0.33 
1.80 ± 0.06 


1.49 ± 0.13 
1.59 ± 0.15 
1.71 ± 0.16 
1.74 ± 0.17 
1.82 ± 0.18 
1.88 ± 0.18 



Table 4: Energies and widths of some prominent lines in the sum spectrum of 
all five detectors determined by our energy calibration and peak fit methods, 
and comparison with the energy given in the literature [84], and the fitted 
dependence of the width as function of energy. 

The background rate in the energy range 2000 - 2080 keV is found to be (0.17 
± 0.01) events/kgy keV {without pulse shape analysis) considering all data as 
background. This is the lowest value ever obtained in such type of experiments. 
The energy resolution at the Q/3/3 value in the sum spectra is extrapolated from 
the strong lines in the spectrum to be (4.00 ± 0.39) keV in the spectra with 
detector 4, and (3.74 ± 0.42) keV (FWHM) in the spectra without detector 4 
(see Fig. 7 and Table 4). The energy calibration of the experiment has an 
uncertainty of 0.20 keV (see Table 4). 
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3 Data Analysis 



We analysed the measured spectra with the following methods: 

1. Bayesian method, which is used widely at present in nuclear and astro- 
physics (see, e.g. [92, 94, 95, 115]). This method is particularly suited for low 
counting rates, where the data follow a Poisson distribution, that cannot be 
approximated by a Gaussian. It has been used also in on earlier double beta 
experiment already [118]. 

2. Feldman-Cousins Method [93, 94]. 

3. Maximum Likelihood Method (see [94, 96]). 

The application of the Bayesian method in our experiment is described in 
[1, 2, 3]. Some numerical examples of the sensitivities of Bayesian and of the 
Maximum Likelihood Methods in the search for events of low statistics are given 
in the next section and in [2, 3, 112, 113]. 

3.1 Numerical Simulations 

To check the methods of analysing the measured data (Bayes and Maximum 
Likelihood Method), and in particular to check the programs we wrote, we have 
generated spectra and lines with a random number generator and performed 
then a Bayes and Maximum Likelihood analysis (see [2, 3]). The length of each 
generated spectrum is 8200 channels, with a line located at bin 5666, the width 
of the line (sigma) being 4 channels (These special values have been choosen 
so that every spectrum is analogue to the measured data). The creation of a 
simulated spectrum is executed in two steps, first the background and second 
the line was created, using random number generators available at CERN (see 
[97] ) . In the first step, a Poisson random number generator was used to calculate 
a random number of counts for each channel, using a mean value of ^=A or 
fi = 0.5 , respectively, in the Poisson distribution 

P{n) = £ e-» (7) 
n\ 

These mean values correspond roughly to our mean background measured in 
the spectra with or without pulse shape analysis. 

In the second step, a Gaussian random number generator was used to calcu- 
late a random channel number for a Gaussian distribution with a mean value of 
5666 (channel) and with a sigma of 4 (channels) . The contents of this channel 
then is increased by one count. This Gaussian distribution filling procedure was 
repeated for n times, n being the number of counts in this line. 

For each choice of fi and n, 100 different spectra were created, and anal- 
ysed subsequently with two different methods: the maximum likelihood method 
(using the program set of [96]) and the Bayes- method. Each method, when 
analysing a spectrum, gives a lower and an upper limit for the number of counts 
in the line for a given confidence level c (e.g. c = 95%) (let us call it confidence 
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area A). A confidence level of 95% means, that in 95% of all cases the true value 
should be included in the calculated confidence area. This should be exactly 
correct when analysing an infinite number of created spectra. When using 100 
spectra, as done here, it should be expected that this number is about the same. 
Now these 100 spectra with a special n and \i are taken to calculate a number 
d, which is the number of that cases, where the true value n is included in the 
resulting confidence area A. 

This number d is given in Table 5 for the results of the two different analysis 
methods and for various values for [i and n. It can be seen, that the Bayes 
method reproduces even the smallest lines properly, while the Maximum Like- 
hood method has some limitations there. 
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Table 5: Results from the analysis of the simulated spectra using a mean back- 
ground of 4 and 0.5 counts, and different line intensities. The number d of cases, 
where the true number of counts in the line (given in the left column), is found 
in the calculated confidence area is given in the second or third column for the 
Bayes method, and in the fourth and fifth column for the maximum likelihood 
method. For details see text. 

Another test has been performed. We generated 1000 simulated spectra 
containing no line. Then the probability has been calculated with the Bayesian 
method that the spectrum does contain a line at a given energy. Table 6 presents 
the results: the first column contains the corresponding confidence limit c (pre- 
cisely the parameter Kg defined earlier), the second column contains the ex- 
pected number of spectra indicating existence of a line with a confidence limit 
above the value c and the third column contain the number of spectra with a 
confidence limit above the value c, found in the simulations. The result un- 
derlines that Ke here is equivalent to the usual confidence level of classical 
statistics. 

We further investigated with the computer-generated spectra the dependence 
of the peak analysis on the width of the energy range of evaluation. Two exam- 
ples are shown in Fig. 8. Here the contents of the simulated peak found with 
the Bayes method is shown as function of the analysis interval given in channels 
(one channel corresponds to 0.36 keV in our measured spectra). The line in the 
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C.L. 


Expected 


Found 


90.0% 


100 ± 31 


96 


95.0% 


50 ± 7 


42 


99.0% 


10 ± 3 


12 


99.9% 


1 ± 1 






Table 6: Number of spectra with a calculated confidence limit above a given 
value. For details see text. 
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Figure 8: Upper part: Analysis of simulated spectrum with Gaussian peak of 
5 events and FWHM of 9.4 channels on a Poisson-distributed background spec- 
trum of 0.5 events/channel, as function of interval of analysis. The middle line 
is the best value. Upper and lower lines correspond to the 68.3% confidence 
limits. Lower part: The same as above, but the peak contains 20 events, the 
background is 4 events/channel. One channel corresponds to 0.36 keV in our 
measured spectra. 

middle is the best-fit value of the method, the upper and lower lines correspond 
to the upper and lower 68.3% confidence limits. In the upper figure the true 
number of counts in the simulated line was 5 events, on a Poisson-distributed 
background of 0.5 events/channel, in the lower figure it was 20 events on a back- 
ground of 4 events/channel. It can be seen, that the analysis gives safely the 
correct number of counts, when choosing an analysis interval of not less than 
40 channels. 

3.2 Analysis of the Full Data 

We first concentrate on the full spectra (see Fig. 4, and Fig. 2 in [1]), without 
any data manipulation (no background subtraction, no pulse shape analysis). 
For the evaluation, we consider the raw data of the detectors. 

The Bayesian peak finding procedure described in [1, 2, 3] determines the 
probability that a line with correct width, and Gaussian shape, exists at given 
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energy E, assuming the rest of the spectra in the interval of analysis as (constant) 
background. It leads for the measured spectra to the result shown on the left 
hand sides of Figs. 10,11. For every energy E of the spectral range 2000 - 
2080 keV, we have determined the probability Kg that there is a line at E. All 
the remainder of the spectrum was considered to be background in this search. 




E(keV) 



Figure 9: Demonstration of the dependence of the relative intensities of the 
weak Bi lines on the location of the impurity in the experimental setup. 
The right side shows two different locations of the impurity (1) relative to the 
detectors. The left side shows the corresponding different relative intensities to 
be expected. Upper part: Source in some distance from the detectors (as in a) 
right part). Lower part: Source in copper cap very near around the detectors 
(as in b) right part). The ordinate shows relative intensities (see [2]). 

The peak detection procedure yields lines at the positions of known [84] weak 
7-lines from the decay of 214 Bi at 2010.7, 2016.7, 2021.8 and 2052.9 keV. The 
lines at 2010.7 and 2052.9 keV are observed at a confidence level of 3.7 and 2.6er, 
respectively. The observed intensities are consistent with the expectations from 
the strong Bi lines in our spectrum, and the branching ratios given in [84] , within 
about the 2a experimental error (see Table 7 and [116]). The expectations here 
are calculated including summing effects of consecutive 7-lines (socalled True 
Coincidence Summing - TCS), by Monte-Carlo simulation of our set-up. Only 
in this way the strong dependence of the relative intensities on the location of 
the impurities in the set-up can be properly taken into account (see Fig. 9). 
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Table 7: 214 Bi is product of the 238 U natural decay chain through j3~ decay of 
214 Pb and a decay of 218 At. It decays to 214 Po by /3~ decay. Shown in this 
Table are the measured intensities of 214 Bi lines in the spectrum shown in Fig.l 
of Ref . [1] in the energy window 2000 - 2060 keV, our calculation of the inten- 
sities expected on the basis of the branching ratios given in Table of Isotopes 
[84] , with and without simulation of the experimental setup, and the intensities 
expected by Aalseth et al. hep-ex/0202018, who do not simulate the setup and 
thus ignore summing of the 7 energies (see [6]). 

*) We have considered for comparison the 3 strongest 214 Bi lines, leaving out 
the line at 1120.287keV (in the measured spectrum this line is partially overim- 
posed on the 1115.55 keV line of 65 Zn). The number of counts in each line have 
been calculated by a maximum-likelihood fit of the line with a gaussian curve 
plus a constant background. 

+) The simulation is performed assuming that the impurity is located in the 
copper part of the detector chamber (best agreement with the intensities of the 
strongest lines in the spectrum). The error of a possible misplacement is not 
included in the calculation. The number of simulated events is 10 8 for each of 
our five detectors. 

**) This result is obtained normalizing the simulated spectrum to the experi- 
mental one using the 3 strong lines listed in column one. Comparison to the 
neighboring column on the right shows that the expected rates for the weak 
lines can change strongly if we take into account the simulation. The reason 
is that the line at 2010.7 keV can be produced by summing of the 1401.50 keV 
(1.55%) and 609.31 keV (44.8%) lines, the one at 2016.7 keV by summing of 
the 1407.98 (2.8%) and 609.31 (44.8%) lines; the other lines at 2021.8 keV and 
2052.94 keV do suffer only very weakly from the summing effect because of the 
different decay schemes. 
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Continues caption of Table 7: H — h) This result is obtained using the number 
of counts for the three strong lines observed in the experimental spectrum and 
the branching ratios from [84] , but without simulation. For each of the strong 
lines the expected number of counts for the weak lines is calculated and then 
an average of the 3 expectations is taken. 

***) Without simulation of the experimental setup. The numbers given here 
are close to those in the neighboring left column, when taking into account that 
Aalseth et al. refer to a spectrum which has only ~11% of the statistics of the 
spectrum shown in Fig.l of Ref. [1]. 



In addition, a line centered at 2039 keV shows up. This is compatible with 
the Q-value [82, 83] of the double beta decay process. We emphasize, that 
at this energy no 7-line is expected according to the compilations in [84] (see 
discussion in section 2). Figs. 10,11 do not show indications for the lines from 
56 Co at 2034. 7keV and 2042 keV discussed earlier [81]. In addition to the line at 
2039 keV we find candidates for lines at energies beyond 2060 keV and around 
2030 keV, which at present cannot be attributed. This is a task of nuclear 
spectroscopy. 

It is of interest, that essentially the same lines as found by the peak scan 
procedure in Figs. 10,11,12, are found [3, 116] when doing the same kind of 
analysis with the best existing natural Ge experiment of D. Caldwell et al. 
[118], which has a by a factor of ~4 better statistics of the background. This 
experiment does, however, not see the line at 2039 keV (see section 3.5). 

Baycsian peak detection (the same is true for Maximum Likelihood peak 
detection) of our data suggests a line at Qpp whether or not one includes detector 
Nr. 4 without muon shield (Figs. 10,11). The line is also suggested in Fig. 12 
after removal of multiple site events (MSE), see below. 




2000 2010 2020 2030 2040 2050 



energy [keV] 




2034 2036 203S 2040 2042 2044 2046 

energy [keV] 



Figure 10: Scan for lines in the full spectrum taken from 1990-2000 with detec- 
tors Nr. 1,2,3,4,5, (Fig. 4), with the Bayesian method of sec. 3.1. The ordinate 
is the probability that a line exists as defined in the text. Left: Energy range 
2000 - 2080 keV. Right: Energy range of analysis ±5cr around Qpp. 



On the left-hand side of Figs. 10,11,12, the background intensity identified 
by the Bayesian procedure is too high because the procedure averages the back- 
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ground over all the spectrum (including lines) except for the line it is trying 
to single out. Inclusion of the known lines into the fit naturally improves the 
background. As example, we show in Fig. 13 the spectrum of Fig. 4 (here in the 
original 0.36 keV binning) with a simultaneous fit of the range 2000 - 2100 keV 
(assuming lines at 2010.78, 2016.70, 2021.60, 2052.94, 2039.0keV and two lines 
at higher energy). The probability for a line in this fit at 2039 keV is 91%. 




energy [keV] energy [keV] 

Figure 1 1 : Left: Probability K e that a line exists at a given energy in the range 
of 2000-2080 keV derived via Bayesian inference from the spectrum taken with 
detectors Nr. 1,2,3,5 over the period August 1990 to May 2000, 46.502 kgy (see 
Fig. 2 from [1].) Right: Result of a Bayesian scan for lines as in the left part of 
this figure, but in the energy range of analysis ± 5<r around Qpp. 




Figure 12: Scan for lines in the single site event spectrum taken from 1995- 
2000 with detectors Nr. 2,3,5, (Fig. 5), with the Bayesian method (as in Figs. 
10,11). Left: Energy range 2000 -2080 keV. Right: Energy range of analysis 
±4.4cr around Q/3/3. 

Finally, on the right-hand side of Figs. 10,11 (and also Fig. 12) the peak 
detection procedure is carried out within an energy interval that seems to not 
contain (according to the left-hand side) lines other than the one at Q/3/3. This 
interval is broad enough (about ± 5 standard deviations of the Gaussian line, 
i.e. as typically used in search for resonances in high-energy physics) for a 
meaningful analysis (see Fig. 8 in section 3.2). We find, with the Bayesian 
method, the probability Kg = 96.5% that there is a line at 2039. OkeV is the 
spectrum shown in Fig. 4. This is a confidence level of 2.1 a in the usual 
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Figure 13: Simultaneous analysis of the spectrum measured with the 76 Ge detec- 
tors Nr. 1,2,3,4,5 over the period August 1990 - May 2000 (54.9813kgy) (same 
as in Fig. 4, but here shown in the 0.36 keV original binning) in the energy 
range 2000 - 2060 keV, with the Maximum Likelihood Method. The probability 
for a line at 2039.0 keV found is this way is 91%. 



language. The number of events is found to be 0.8 to 32.9 (7.6 to 25.2) with 
95% (68%) c.L, with best value of 16.2 events. For the spectrum shown in Fig. 
2 in [1], we find a probability for a line at 2039.0 keV of 97.4% (2.2 a). In this 
case the number of events is found to be 1.2 to 29.4 with 95% c.L. It is 7.3 to 
22.6 events with 68.3% c.L. The most probable number of events (best value) 
is 14.8 events. These values are stable against small variations of the interval 
of analysis, as expected from Fig. 8 in section 3.2. For example, changing the 
lower and upper limits of the interval of analysis between 2030 and 2032 and 
2046 and 2050 yields consistently values of between 95.3 and 98.5% (average 
97.2%) for the spectrum of Fig.2 of [1]. 

We also applied the Feldman-Cousins method recommended by the Particle 
Data Group [93, 94]. This method (which does not use the information that the 
line is Gaussian) finds a line at 2039 keV on a confidence level of 3.1 a (99.8% 
c.L). 



3.3 Analysis of Single Site Events Data 

From the analysis of the single site events (Fig. 5), we find after 28.053 kgy of 
measurement 9 SSE events in the region 2034.1 - 2044.9 keV (± 3tr around Qpp) 
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(Fig. 14). Analysis of the single site event spectrum (Fig. 5), as described in 
section 3.1, shows again evidence for a line at the energy of Qpp (Fig. 12). An- 
alyzing the range of 2032 - 2046 keV, we find the probability of 96.8% that there 
is a line at 2039.0 keV. We thus see a signal of single site events, as expected for 
neutrinoless double beta decay, at the value obtained in the precision ex- 
periment of [82]. The analysis of the line at 2039.0 keV before correction for the 
efficiency yields 4.6 events (best value) or (0.3 - 8.0) events within 95% c.l. ((2.1 
- 6.8) events within 68.3% c.l.). Corrected for the efficiency to identify an SSE 
signal by successive application of all three PSA methods, which is 0.55 ± 0.10, 
we obtain a 0^/3/3 signal with 92.4% c.l.. The signal is (3.6 - 12.5) events with 
68.3% c.l. (best value 8.3 events). Thus, with proper normalization concerning 
the running times (kgy) of the full and the SSE spectra, we see that almost the 
full signal remains after the single site cut (best value), while the 214 Bi lines 
(best values) are considerably reduced (see Fig. 15). 

' l Range of Bi lines + 0v(S(S Signal Calibration Lines 




2010 2016 2021 2053 2039 609.6 1763.9 1592 2103 2614 

' ^ > ' ' " ' 



Figure 15: Relative suppression ratios: Remaining intensity after pulse shape 
analysis compared to the intensity in the full spectrum. Right: Result of a cali- 
bration measurement with a Th source - ratio of the intensities of the 1592 keV 
line (double escape peak, known to be 100% SSE), set to 1. The intensities of 
the 2203 keV line (single escape peak, known to be 100% MSE) as well as of 
the 2614 keV line (full energy peak, known to be dominantly MSE) are strongly 
reduced (error bars are Ha). The same order of reduction is found for the 
strong Bi lines occurring in our spectrum (see Fig. 3) - shown in this figure are 
the lines at 609.4 and 1763.9 keV. Left: The lines in the range of weak statistics 
around the line at 2039 keV (shown are ratios of best fit values). The Bi lines 
are reduced compared to the line at 2039 keV (set to 1), as to the 1592 keV SSE 
Th line. 

We have used a 23& Th source to test the PSA method. We find the reduction 
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of the 2103 keV and 2614 keV 228 Th lines (known to be multiple site or mainly 
multiple site), relative to the 1592 keV 228 Th line (known to be single site), 
shown in Fig. 15. This proves that the PSA method works efficiently. Essen- 
tially the same reduction as for the Th lines at 2103 and 2614 keV and for the 
weak Bi lines is found for the strong 214 Bi lines (e.g. at 609.6 and 1763.9 keV 
(Fig. 15)). 

The Feldman-Cousins method gives a signal at 2039.0 keV of 2.8 a (99.4%). 
The possibility, that the single site signal at 2039.0 keV is the double escape 
line corresponding to a (much more intense!) full energy peak of a 7- line at 
2039+1022=3061keV is excluded from the high-energy part of our spectrum. 

3.4 Comparison with Earlier Results 

We applied the same methods of peak search as used in sections 3.3, 3.4, to the 
spectrum measured in the Ge experiment by Caldwell et al. [118] more than a 
decade ago. These authors had the most sensitive experiment using natural Ge 
detectors (7.8% abundance of 76 Ge). With their background being a factor of 
9 higher than in the present experiment, and their measuring time of 22.6 kg y, 
they have a statistics for the background larger by a factor of almost 4 in their 
(very similar) experiment. This allows helpful conclusions about the nature of 
the background. 

The peak scanning finds [3, 117] (Fig. 16) indications for peaks essentially 
at the same energies as in Figs. 10,11,12. This shows that these peaks are not 
fluctuations. In particular it sees the 2010.78, 2016.7, 2021.6 and 2052.94keV 
214 Bi lines. It finds, however, no line at Qpp (see also Fig. 17). This is consistent 
with the expectation from the rate found from the HEIDELBERG-MOSCOW 
experiment. About 16 observed events in the latter correspond to 0.6 expected 
events in the Caldwell experiment, because of the use of non-enriched material 
and the shorter measuring time. 

Another Ge experiment (IGEX) using 9 kg of enriched 76 Ge, but collecting 
since beginning of the experiment in the early nineties till shutdown in end of 
1999 only 8.8 kgy of statistics [125], because of this low statistics also naturally 
cannot see any signal at 2039 keV. It sees, however, also the line around 2030 keV 
indicated in Figs. 10,11,12 from our data. 

3.5 Proofs and Disproofs 

The result described in sections 3.3-3.5 has been questioned in some papers 
(Aalseth et al, hep-ex/0202018; Feruglio et al., Nucl. Phys. B 637(2002)345; 
Zdesenko et al., Phys. Lett. B 546(2002) 206, and Kirpichnikov, talk at Meeting 
of Physical Section of Russian Academy of Sciences, Moscow, December 2, 2002, 
and priv. communication, Dec. 3, 2002.) These claims against our results are 
incorrect in various ways. 



24 



S.ioo 
















■i so 

al 




(1 
















GO 


























40 


























20 




, , 1 , 


, , 1 , 


r 




r 















2000 2010 2020 2030 2040 2050 2060 2070 2080 2090 21 OO 

etiergy(KeV) 




2000 2010 2020 2030 2040 2050 2060 2070 2080 2090 21 OO 

etiergy(KeV) 

Figure 16: Peak scanning of the spectrum measured by Caldwell et al. [118] 
with the Maximum Likelihood method (upper part), and with the Bayesian 
method (lower part) (as in Fig. 10,11,12) (see [116]). 



The arguments in the first two of these papers can be easily rejected. We have 
published this in hep-ph/0205228. In particular their estimates of the intensities 
of the 21A Bi lines are simply wrong, because the summing effect of the energies 
of consecutive gamma lines (known as True Coincidence Summing - TCS in 
nuclear spectroscopy) was not taken into account (see Table 7 in the present 
paper, and in [3]). Further none of the papers and in has performed a Monte 
Carlo simulation of our setup, which is the only way to come to quantitative 
statements. 

The paper by Zdesenko et al. starts from an arbitrary assumption, namely 
that there are lines in the spectrum at best only at 2010 and 2053 keV. This 
contradicts to the experimental result, according to which there are further 
lines in the spectrum (see Figs. 10,11,12,16 in the present paper). In this way 
they artificially increase the background in their analysis and come to wrong 
conclusions. 

All three of these papers, when discussing the choice of the width of the 
search window, ignore the results of the statistical simulations we give in hep- 
ph/0205228 and we have published in [2, 3] (see also Fig. 8 in the present 
paper). 
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Figure 17: Analysis of the spectrum measured by D. Caldwell et al. [118], with 
the Maximum Likelihood Method, in the energy range 2000 - 2060 keV (as in 
Fig. 13) assuming lines at 2010.78, 2016.70, 2021.60, 2052.94, 2039.0 keV. In 
contrast to Fig. 13 no indication for a signal at 2039 kev is observed (see text). 



Kirpichnikov states that from his analysis he clearly sees the 2039 keV line 
in the full (not pulse shape discriminated) spectrum with > 99% c.l. He claims 
that he does not see the signal in the pulse shape spectrum. One reason to see 
less intensity there certainly is that in this case he averages for determination 
of the background over the full energy range without allowing for any lines. His 
result is in contradiction with the result we obtain under the same assumption 
of assuming just one line (at Q/3/3) and a continuous background (see Fig. 12 of 
this paper). 

4 Half-life of the Neutrinoless Mode and 
Effective Neutrino Mass 

Having shown that the signal at Q/3/3 consists of single site events and is not 
a 7 line, we translate the observed number of events into half-lives for the 
neutrinoless double beta decay. We give in Table 8 conservatively the values 
obtained with the Bayesian method and not those obtained with the Feldman- 
Cousins method. 

We derive from the data taken with 46.502 kg y the half-life T°J 2 = (0.8 - 
18.3) x 10 25 y (95% c.l.). The analysis of the other data sets, shown in Table 8 
confirm this result. Of particular importance is that we see the 0^/3/3 signal in 
the single site spectrum. 

The result obtained is consistent with the limits given by all other double 
beta experiments - which still reach less sensitivity (see Figs. 18,19). The most 
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Figure 18: Present evidence for Ou/3/3 decay from data of the HEIDELBERG- 
MOSCOW experiment and the potential of present and future (3(3 experiments. 
Vertical axis - half life limits (only HEIDELBERG-MOSCOW gives a value) in 
years, horizontal - isotopes used in the various experiments. 



sensitive experiments following the HEIDELBERG-MOSCOW experiment are 
the geochemical 128 Te experiment with T^ 2 > 2(7.7) x 10 24 y (68% c.L), [120] 
the 136 Xe experiment by the DAMA group with T°^ 2 > 1.2 x 10 24 y (90% 
c.l.)[122], a second 76 Ge experiment with T% > 1.2 x 10 24 y [119] and a nat Ge 

experiment with > 1 x 10 24 y [118, 119]. Other experiments are already 
about a factor of 100 less sensitive concerning the 0v/3(3 half-life: the Gotthard 
TPC experiment with 136 Ae yields [123] T°% > 4.4 x 10 23 y (90% c.L) and the 
Milano Mibeta cryodetector experiment [128] T°^ 2 > 1.44 x 10 23 y (90% c.L). 

Another experiment [125] with enriched 76 Ge, which has stopped operation 
in 1999 after reaching a significance of 8.8 kg y, yields (if one believes their 
method of 'visual inspection' in their data analysis), in a conservative analysis, 
a limit of about T°^ 2 > 5 x 10 24 y (90% c.L). The 128 Te geochemical experiment 
yields (m„) < 1.1 eV (68 % c.L) [120], the DAMA 136 Ae experiment (m v ) < 
(1.1 — 2.9) eV [122] and the 130 Te cryogenic experiment yields (m v ) < 1.8 eV 
[128]. 

Concluding we obtain, with more than 95% probability, first evidence for 
the neutrinoless double beta decay mode. 

As a consequence, at this confidence level, lepton number is not conserved. 
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Figure 19: Present evidence for Ou/3/3 decay from data of the HEIDELBERG- 
MOSCOW experiment and the potential of present and future (3(5 experiments. 
Vertical axis - effective neutrino mass limits (only HEIDELBERG-MOSCOW 
gives a value) in eV, horizontal - isotopes used in the various experiments. 



Further our result implies that the neutrino is a Majorana particle. Both of 
these conclusions are independent of any discussion of nuclear matrix elements. 

The matrix element enters when we derive - under the assumption that the 
Qv(3(3 amplitude is dominated by the neutrino mass mechanism, a value for the 
effective neutrino mass. If using the nuclear matrix element from [39, 38], we 
conclude from the various analyses given above the effective mass (to) to be 
(to) = (0.11 - 0.56) eV (95% c.L), with best value of 0.39eV (see Tabl. 8). 
Allowing conservatively for an uncertainty of the nuclear matrix elements of ± 
50% (for detailed discussions of the status of nuclear matrix elements we refer 
to [27, 28, 40, 43, 8, 47, 45, 46]) this range may widen to (to) = (0.05 - 0.84) eV 
(95% c.L). In Table 9 we demonstrate the situation of nuclear matrix elements 
by showing the neutrino masses deduced from different calculations (see [3] ) . It 
should be noted that the value obtained in Large Scale Shell Model Calculations 
[41] is understood to be too large by almost a factor of 2 because of the two 
small configuration space, (see, e.g. [43]), and that the second highest value 
given (from [42]), has now been reduced to 0.53 eV [47]. The recent studies 
by [45, 46] yield an effective mass of (0.44 - 0.52) eV. We see that the early 
calculations [39] done in 1989 agree within less than 25% with the most recent 
values. 
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Significan- 
ce [kgy] 


Detectors 


T?/ 2 y 


(to) eV 


Conf. 
level 


54.9813 


1,2,3,4,5 


(0.80- 35.07) x 10 25 
(1.04- 3.46) x 10 25 
1.61 x 10 25 


(0.08 - 0.54) 
(0.26 - 0.47) 
0.38 


95% c.l. 
68% c.l. 
Best Value 


46.502 


1,2,3,5 


(0.75 - 18.33) x W Zb 
(0.98 - 3.05) x 10 25 
1.50 x 10 25 


(0.11 - 0.56) 
(0.28 - 0.49) 
0.39 


95% c.l. 
68% c.l. 
Best Value 


28.053 


2,3,5 SSE 


(0.88 - 22.38) x 10 25 
(1.07- 3.69) x 10 25 
1.61 x 10 25 


(0.10 - 0.51) 
(0.25 - 0.47) 
0.38 


90% c.l. 
68% c.l. 
Best Value 



Table 8: Half- life for the neutrinoless decay mode and deduced effective neutrino 
mass from the HEIDELBERG-MOSCOW experiment. 



M 0u from 


[38, 39] 


[40] 


[37] 


[41] 


[42] 


[44] 


[47] 


[45, 46] 


(to) cV 95% 


0.39 


0.37 


0.34 


1.06 


0.87 


0.60 


0.53 


0.44 - 0.52 



Table 9: The effect of nuclear matrix elements on the deduced effective neutrino 
masses. Shown are the neutrino masses deduced from the best value of Tjj 2 = 
1.5 x 10 25 y determined in this work with 97% c.l. when using matrix elements 
from various calculations and a phase factor of F® v — 6.31 x 10~ 15 

In the above conclusion for the effective neutrino mass, it is assumed that 
contributions to 0vf3(3 decay from processes other than the exchange of a Ma- 
jorana neutrino (see, e.g. [8, 98] and section 1) arc negligible. It has been 
discussed, however, recently [61] that the possibility that 0v/3(3 decay is caused 
by R-parity violation, may experimentally not be excluded, although this would 
require making R-parity violating couplings generation-dependent. 

Assuming other mechanisms to dominate the Qv(3f3 decay amplitude, the re- 
sult allows to set stringent limits on parameters of SUSY models, leptoquarks, 
compositeness, masses of heavy neutrinos, the right-handed W boson and pos- 
sible violation of Lorentz invariance and equivalence principle in the neutrino 
sector. For a discussion and for references we refer to [8, 30, 23, 33, 34]. 
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Figure 20: The impact of the evidence obtained for neutrinoless double beta 
decay in this paper (best value of the effective neutrino mass (m) = 0.39 eV, 
95% confidence range (0.05 - 0.84) eV - allowing already for an uncertainty of the 
nuclear matrix element of a factor of ± 50%) on possible neutrino mass schemes. 
The bars denote allowed ranges of (m) in different neutrino mass scenarios, still 
allowed by neutrino oscillation experiments. Hierarchical models are excluded 
by the new Ov(3(3 decay result. Also shown are the expected sensitivities for the 
future potential double beta experiments CUORE, MOON, EXO and the 1 ton 
and 10 ton project of GENIUS [100] (from [5]). 

5 Consequences 

With the value deduced for the effective neutrino mass, the HEIDELBERG- 
MOSCOW experiment excludes several of the neutrino mass scenarios allowed 
from present neutrino oscillation experiments (see Fig. 20) - allowing only for 
degenerate mass scenarios, and an inverse hierarchy 3v and Av- scenario. For 
details we refer to [5]. In particular, hierarchical mass schemes are excluded at 
the present level of accuracy. 

According to [55] a global analysis of all solar neutrino data including the 
recent SNO neutral-current rate selects the Large Mixing Angle (LMA) at the 
90% c.L, however, the LOW solution is also viable, with 0.89 goodness of fit. 
Recently the LMA solution has been supported by evidence for reactor antineu- 
trino disappearance [56]. 

Assuming the degenerate scenario to be realized in nature we fix - according 
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to the formulae derived in [7] - the common mass eigenvalue of the degenerate 
neutrinos to m = (0.05 - 3.4) cV. Part of the upper range is already excluded 
by tritium experiments, which give a limit of m < 2.2 eV, or 2.8 eV (95% c.l.) 
[99]. The full range can only partly (down to ~ 0.5 eV) be checked by future 
tritium decay experiments, but could be checked by some future (3(3 experiments 
(see, e.g. [31, 8, 100]). The deduced 95% interval for the sum of the degenerate 
neutrino masses is consistent with the range for Sl v deduced from recent cosmic 
microwave background measurements and large scale structure (rcdshift) sur- 
veys, which still allow for a J^i m i 5= 4.4 eV [109, 108]. The range of (m) fixed 
in this work is, already now, in the range to be explored by the satellite experi- 
ments MAP and PLANCK [107, 110] (see Fig. 21). It lies in a range of interest 
for Z-burst models recently discussed as explanation for super-high energy cos- 
mic ray events beyond the GKZ-cutoff [103, 104, 72]. Finally, the deduced best 
value for the mass is consistent with expectations from experimental fi — > e-f 
branching limits in models assuming the generating mechanism for the neutrino 
mass to be also responsible for the recent indication for an anomalous magnetic 
moment of the muon [101]. A recent model with underlying A4 symmetry for 
the neutrino mixing matrix (and the quark mixing matrix) also leads to de- 
generate neutrino masses consistent with the present experiment [102]. This 
model succeeds to consistently describe the large (small) mixing in the neutrino 
(quark) sector. 

The neutrino mass deduced leads to 0.002 < Q u h 2 < 0.1, and thus may allow 
neutrinos to still play an important role as hot dark matter in the Universe (see 
also [60, 35]). 

6 Future of (3(3 Experiments - GENIUS and 
Other Proposals 

With the HEIDELBERG-MOSCOW experiment, the era of the small smart 
experiments is over. New approaches and considerably enlarged experiments 
(as discussed, e.g. in [8, 31, 32, 100, 105]) will be required in future to fix the 
neutrino mass with higher accuracy. 

Since it was realized in the HEIDELBERG-MOSCOW experiment, that the 
remaining small background is coming from the material close to the detec- 
tor (holder, copper cap, ...), elimination of any material close to the detector 
will be decisive. Experiments which do not take this into account, like, e.g. 
CUORE [128, 129], and MAJORANA [126], will allow only rather limited steps 
in sensitivity. Further there is the problem in cryodetectors that they cannot 
differentiate between a [3 and 7 signal, as this is possible in Ge experiments. 

Another crucial point is - see eq. (6) - the energy resolution, which can 
be optimized only in experiments using Germanium detectors or bolometers. 
It will be difficult to probe evidence for this rare decay mode in experiments, 
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Figure 21: Double beta decay observable (m) and oscillation parameters: The 
case for degenerate neutrinos. Plotted on the axes are the overall scale of neu- 
trino masses mo and mixing tan 2 8 12 . Also shown is a cosmological bound 
deduced from a fit of CMB and large scale structure [107] and the expected 
sensitivity of the satellite experiments MAP and PLANCK. The present limit 
from tritium (3 decay of 2.2 eV [99] would lie near the top of the figure. The 
range of (to) fixed by the HEIDELBERG-MOSCOW experiment is, in the case 
of small solar neutrino mixing, already in the range to be explored by MAP and 
PLANCK. 

which have to work - as result of their limited resolution - with energy windows 
around of up to several hundreds of keV, such as NEMO III [130], EXO 
[127], CAMEO [124]. 

Another important point is (see eq. 6), the efficiency of a detector for detec- 
tion of a (3(3 signal. For example, with 14% efficiency, a potential future 100 kg 
82 Se NEMO experiment would be equivalent only to a 10 kg experiment (not 
talking about the energy resolution). 

In the first proposal for a third generation double beta experiment, the GE- 
NIUS proposal [31, 23, 100], the idea is to use 'naked' Germanium detectors 
in a huge tank of liquid nitrogen. It seems to be at present the only proposal, 
which can fulfill both requirements mentioned namely to increase the detector 
mass and to simultaneously reduce the background drastically. The potential 
of GENIUS is together with that of some later proposals indicated in Fig. 20. 
GENIUS would - with only 100 kg of enriched 76 Ge - increase the confidence 
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level of the present pulse shape discriminated signal to 4cr within one year, and 
to 7a within 3 years of measurement (a confirmation on a 4a level by the MA JO- 
RANA project would need ~ 230 years, while the CUORE project would need 
(ignoring the problem of identification of the signal as (3(3 signal) 3700 years). 
With ten tons of 76 Ge GENIUS should be capable to investigate also, whether 
the neutrino mass mechanism or another mechanism (see, e.g. [8]) is dominating 
the decay amplitude. A GENIUS Test Facility is at present under construction 
in the GRAN SASSO Underground Laboratory [132, 134]. 

7 Conclusion 

The status of present double beta decay search has been presented, and re- 
cent evidence for a non-vanishing Majorana neutrino mass obtained by the 
HEIDELBERG-MOSCOW experiment about one year ago has been described 
and outlined. Future projects to improve the present accuracy of the effective 
neutrino mass have been briefly discussed. The most sensitive of them and 
perhaps at the same time most realistic one, is the GENIUS project. GENIUS 
is the only of the new projects which simultaneously has a huge potential for 
cold dark matter search, and for real-time detection of low-energy neutrinos (see 
[30, 36, 131, 33, 133, 106, 34]). 
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